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ABSTRACT 

The properties of gas in the halos of galaxies constrain global models of the interstellar 
medium. Kinematical information is of particular interest since it is a clue to the origin 
of the gas. Until now mostly massive galaxies have been investigated for their halo 
properties. 

Here we report on deep H i and Ha observations of the edge-on dwarf galaxy UGC 
1281 in order to determine the existence of extra-planar gas and the kinematics of 
this galaxy. This is the first time a dwarf galaxy is investigated for its gaseous halo 
characteristics. We have obtained Ha integral field spectroscopy using PPAK at Calar 
Alto and deep H i observations with the WSRT of this edge-on dwarf galaxy. These 
observations are compared to 3D models in order to determine the distribution of H i 
in the galaxy. 

We find that UGC 1281 has Ha emission up to 25"(655 pc) in projection above 
the plane and in general a low Ha flux. Compared to other dwarf galaxies UGC 1281 
is a normal dwarf galaxy with a slowly rising rotation curve that flattens off at 60 
km and a central depression in its Hi distribution. Its Hi extends 70"(1.8 kpc) 
in projection from the plane. This gas can be explained by either a warp partially in 
the line-of-sight or a purely edge-on warp with rotational velocities that decline with 
a vertical gradient of 10.6 ± 3.7 km s~^ kpc~^. The line-of-sight warp model is the 
preferred model as it is conceptually simpler. In either model the warp starts well 
within the optical radius. 

Key words: UGC 1281, Galaxies: halos, Galaxies: ISM, Galaxies: spiral, Galaxies: 
structure, Galaxies: kinematics and dynamics. Galaxies: dwarf 



1 INTRODUCTION 

The discovery of a very extended H I halo in NGC 891 
has shown that 21 cm observations are a powerful tool 
to study the structu re and kinematics of gaseous halos. 
Since this discovery jSwaters et al.l Il997 ) many investiga - 



tions have followed studying this ( Oosterloo et al.l 120071 ') 
and other galaxies. Investigations at 21 cm, as well as 
other wavelengths, have shown that NGC 891 is not a 



special case in having a gaseous halo JSchaao et al.l 


200C 


Lee et al.ll200ll:lRossa & Dettmaill2003l: 


Barbieri et al.l 


2005 


Westmeier et al.l 20051: Boomsma et al.l 


20081). 



It is thought that a large part of this h alo gas is brought 
up from the disk b y galactic fountains l|Shapiro fc Field! 
1 19761 : iBregmanI Il980l ) . In this model supernovae explosions 
expel the hot hydrogen gas from the disk into the halo. In 
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the halo the gas cools down and condenses before returning 
to the disk. Studies of massive galaxies have shown that this 
mechanism is very likely at work in several galaxies, either 
inferred from the fact t hat the extra-planar H I is concen- 
trated towards the disk jOosterloo et al.ll2007l ) or that the 
majority o f high velocity clouds are located near the bright 
inner disk (|Boomsma et al.ll2008l ). 

The gas in the halo of NGC 891 shows rotational ve- 
locities that are lower then the v elocities of the g as in the 
disk. This 'lag ' is observed in H I ('Fraternali et al."2005l) as 
well as in Ha (|Heald et al.ll2006l : iKamphuis ct al. 2007) and 
determined to be linear with a magnitude of ^-^16 km s~^ 
kpc~^ in the vertical direction. 

Simulations have shown that the lag cannot be expl ained 
by galactic fountains alone (|Fraternali fc Binnevll200g ). To- 
gether with other observational facts, such as the low metal- 
licity of the High Velocity Clouds (HVCs) and the presence 
of filaments and other irregular gaseous structures in the 
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halo l|Sancisi et al.|[2008l ). this has given rise to the idea 
that the halo must also be partially formed from gas that has 
been accreted from the Intergalactic Medium (IGM) or com- 
panion galaxies. However, the gradients by the mselves do 
not necessarily require ac cretion from the IGM l|Beniaminl 
l2002l : lBarnab6i et aLllioO^ l. 

Even though this lagging b ehavior is ob s erved in other 
galaxies then NGC 891 -e.g. i Heald et all ((20071) studied 
two other large galaxies and found a lag in both of them- 
it is not known whether it occurs in all gaseous halos. 
Until now, kinematical studies have focussed on massive 
galaxies and one Low S urface Brightness (LSB) galaxy 
l|Matthews fc WoodllioO^ I but dwarf galaxies have not yet 
been investigated for a kinematic lag. 

Here we present deep H l and Ha observations of the 
edge-on dwarf galaxy UGC 1281 and analyze the kinemat- 
ics of the H I and Ha in and above the plane. The 21 cm line 
emission was observed with the Westerbork Synthesis Radio 
Telescope (WSRT) and the Ha with the Integral Field Unit 
(IFU) PPAK on the 3.5 m telescope at Calar Alto. 

The edge-on orientation of UGC 1281 provides us with 
an excellent opportunity to study the vertical structure of 
gas in a dwarf galaxy. Our H I observations of UGC 1281 are 
some of the deepest observations ever conducted on a dwarf 
galaxy, and never before has an edge-on dwarf galaxy been 
observed to this depth. 

UGC 1281 is a nearby e dge-on dwarf galaxy (Mb=- 
15.8. ISwaters fc Balcelld (|2002l )) with a systemic velocity of 
156 k m s~^ at a distance of 5.4 Mpc l|Karachentsev et al.l 
bOO^I. The galaxy ha s a angular diameter of 4.46' (D25) 
( de Vaucouleurs et al.l ll992^) on the s ky which would trans- 
late to 7 kpc at the given distance. Ivan Zeel (|200Cll . I2OOII ) 
measured the star formation rate in UGC 1281 to be very 
low (0.0084 M0 yr.~^). However, this should be consid- 
ered as a lower limit since the Ha flux was not corrected 
for internal extinction. These observati ons are supported by 
the non-detections in radio continuum jHummel et al.iri99ll ) 
and IRAS. Indeed. iRossa fc Dettm^l \2Q0± observed UGC 
1281 as a part of their extra-planar diffuse ionized gas (DIG) 
survey and could not detect any extra-planar Ha in this 
dwarf galaxy. All these observations indicate a very low star 
formation rate in UGC 1281. 

Under the assumption that gaseous halos are created 
by processes related to star formation, little to no halo 
is expected when the SFR is low. However, in a dwarf 
galaxy gas might easily escape from the disk due to the 
weak potential. Also, even though small galaxies are not 
expected to accrete bary onic matter at lower redshifts 
(|Hoeft fc GottloebeJlioiol ). UGC 1281 is in the transition 
region, between accreting and non-accreting galaxies, and 
therefore a modest amount of lagging extra-planar gas might 
be present if the vertical gradient is predominantly formed 
by the accretion of matter. 

This article is structured as followed. In § [2] we will de- 
scribe the observations and data reduction. §[3] will contain 
the results of the observations followed by a presentation of 
our models in § [4] We will discuss our results in § [5] and 
summarize in §(6] 



Parameter Value 



Observation date 2004 September 

Total length of observation (hr) 4 X 12 

Velocity center of band (km s^^) 156 

Total bandwidth (MHz) 10 

Channels in obs. 1024 

Channel sep. in obs. (kHz) 9.77 

Channels in final cube 61 



Vel. res. after Hanning smoothing (km s ^) 8.24 



Table 1. Log of the Hi observations. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Radio Data 

The 21 cm line emission, or Hi, observations were obtained 
with the WSRT during four nights in September 2004. In 
total 4 complete 12 hr observations were performed using 
the Maxi-Short configuration. This configuration gives the 
optimum performance for imaging extended sources. This 
is because it provides a good sampling of the inner U — V 
plane with a shortest base line of 36 m. The longest baseline 
is 2754 m. An overview of the observational parameters is 
given in Table [T] 

The data we re reduced using the MIRIAD package 
ISauh et al.l 1 19951 ). Care was taken not to include data 
affected by solar or other interference. Before and after each 
12 hr observation a calibration source (3C147 and CTD93) 
was observed, thus enabling us to determine the spectral 
response of the telescope. During each 12 hr observation 
no additional (phase) calibration sources were observed, 
as is standard practice with the WSRT. Due to the large 
bandwidth, the time variations can be determined from 
self-calibration of detected sources in the continuum image 
constructed from the channels free of line emission. In 
this case UGC 1281 itself was not detected in continuum 
emission. 

After the reduction and calibration additional 
analysis was performed with the GIPSY package 
(Ivan der Hulst et"allll992l ). The final cube was reduced to 
61 velocity channels with a velocity spacing of 4.12 km 
s^^, which results in a velocity resolution after Hanning 
smoothing of 8.24 km s ~^ . The original cube has a spatial 
resolution of 25.0" x 13.4". This cube was smoothed to two 
cubes with circular beams of 26"(0.68 kpc) and 60"(1.57 
kpc) to avoid beam orientation effects and to detect low 
level emission (see Table [Sjl. 

For the calculation of the minimum detectable column 
in Table[2]we used a velocity width of 16.5 km s~^, because 
real emission will never appear in a single channel. 
The position of the cente r of the galaxy, as determined by 
ISwaters fc Balcelli (|2002l ') from _R-band photometry, was 
set to zero i n the three cubes and they were rotated by 50 ° 
(PA = 40° . ISwaters fc Balcellj \2mm to orient the major 
axis of the galaxy parallel to the x-axis of the image. 
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Parameter 


Full Resolution 


Circular Beam 


Low Resolution 


Spatial resolution (") 


25.0 X 13.4 


26.0 X 26.0 


60.0 X 60.0 


Beam size (kpc) 


0.65 X 0.35 


0.68 X 0.68 


1.57 xl.57 


rms noise per channel (mjy beam~^) 


0.44 


0.50 


0.72 


Minimum detectable column density (3cr; cm~^) 


7.1 X 10i3 


4.0 X 10" 


1.1 X 10" 



Table 2. Parameters of the Hi Data Cubes. 



Parameter Field 1&2 Field 3 



Name North&Center South 

Observation date 2006 September same 

Exposure time (s) 3600 1800 

Central wavelength (A) 6273 6273 

Total bandwidth (A) 1628 1628 

Channels in obs. 1050 1050 

Channel sep. in obs. (A) 1.55 1.55 

Channel separation (km s~^) at Ha 70.8 70.8 



Table 3. Log of the Ha observations. 



in the three fields were ahgned properly, this procedure re- 
sulted in a position error less than 1", which is smaller than 
the 2.7" fiber size. 

Our reduced and calibrated data were further ana- 
lyzed with a combination of IDL programs prod u ced b y 
the SAURON collaboration (e.g. Ide Zeeuw (|2002l ')') 
and ourselves. The data were Vo ronoi binned with 
the y oronoi_2d_binning IDL program (jCaDDellari fc Copinl 
l2003h to obtain a higher S/N ratio (S/N > 20). Like the H I, 
the central position of the galaxy was set to (0,0) and the 
galaxy was rotated by 50° to align the major axis with the 
X-axis. 



2.2 IFU data 

The Hq was observed with the PP AK integral field unit on 
the 3.5 m telescope at Calar Alto (|Kelz et allbood ). Three 
positions were observed in one night (see Table [3]). At each 
position several exposures were taken. Two positions were 
observed for 3 x 1200 s (North and Center) and one for 3 x 
600 s (South). The lower exposure time of the Southern field 
was due to twilight. The pointings of the North and South 
fields were shifted by ~ ± 76", compared to the Center, 
along an axis with a PA slightly offset (~ 10° offset) from 
the PA of the galaxy. This was done to ensure a good cov- 
erage of extra-planar Ha. The spectra cover a wavelength 
range ~5500 - 7000 A with a resolution of 4.1 A(187 km 
at Ha). The conditions were partially photometric with 
a mean seeing of 1.3", which is much less than the fiber size 
(2.7"). 

For the reduction of the PPAK spectra th e IRAF pack - 
age was used. The steps in the dohydra Guide (|ValdeJl995l l 
were followed manually to ensure complete control over the 
reduction. The only deviation taken from the steps as de- 
scribed in this guide was that we traced the apertures on 
the science frames themselves. This was done because there 
was more than enough continuum emission in each spectrum 
and there was no need to introduce additional errors due to 
shifting the apertures. 

For the initial wavelength calibration 7 lines of a 
HeNe+ThAr lamp were used. Afterwards a fine tuning cali- 
bration was performed with 5 sky lines in each science frame. 
For the sky line subtraction a total of 95 sky fibers from the 
three pointings were used. All sky fibers were checked for 
inconsistencies and, except for sky fibers with the galaxy in 
their field of view, none were found. 

After the reduction and wavelength calibration the dif- 
ferent exposures of each field were combined by calculating 
the median of the separate exposures. No clipping was ap- 
plied. The three fields were then positioned into one field by 
overlaying the continuum images of the fields on top of a R- 
band DSS image. By taking special care that the five stars 



3 RESULTS 

Here, we first discuss the distribution and kinematics of the 
hydrogen that can be directly derived from the data. In Sec- 
tion [S] we will compare the distributions and kinematics of 
the ionized and neutral gas, highlight the differences and 
similarities, and discuss possible interpretations of these dis- 
tributions. 



3.1 Gas distribution 

3.1.1 Ha Distribution 

Fig. [1] shows the Ha fiux obtained with the PPAK IFU 
instrument. The flux in each bin was determined by fltting 
a Gaussian, with the IDL routine GAUSSFIT, to an average 
of the spectra in each bin. The area under this Gaussian is 
then considered to be the fiux in each bin. All the fitted 
Gaussians were inspected by eye, and none of the fitted 
lines showed significant deviations from Gaussian. For 
display purposes the log values of the fiux are displayed. 
For the same reason the bins without emission are simply 
taken away from the data and thus not displayed. Overlaid 
on the image are the contours of an integrated H I (zeroth 
moment) map. 

The first thing we see from Fig. [1] is that the ionized 
hydrogen is mostly located in 5 or 6 distinct peaks but 
that there is low level emission almost everywhere in the 
central field of view. The distinct peaks can most likely be 
associated with large H 1 1 regions in the galaxy whereas the 
low level emission is indicating a diffuse ionized component. 
The two innermost peaks at the NE and SW side of the 
center of the galaxy might be caused by a ring-like structure 
or central depression in the distribution. This seems to be 
implied by their equal distance from the center; however, 
they could also be normal H 1 1 regions. 

A warp is clearly seen in the H 1 1 regions as well as the 
diffuse emission. The peak of the emission starts to deviate 
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Figure 1. Integrated Ha flux map of UGC 1281. The map was constructed by taking the log of the area under the fitted Gaussian in 
each spectrum. The red lines indicate the full field of view of the PPAK observations. Overlaid in black are the contours of an integrated 
Hi (zeroth moment) map. The contours are at 1.5, 3, 6, 12, 24, 48, 96, 192, 273 a levels of the data where a = 2.2 mjy beam"^. The 
arrows indicate the positions of the velocity cuts parallel to the minor axis in Fig. |6] 



from the major axis at a radial distance of ^ 50"(1.3 kpc). 
In the South West (SW) of the galaxy we see an exception 
to this behavior with a large ionized hydrogen peak below 
the major axis, where the diffuse emission still seems to 
be mostly on the major axis. This seems to indicate an 
H 1 1 region that is either somewhat offset from the plane 
or located in the outskirts of the galaxy. If the warp is 
partially along the line of sight the outer parts will not only 
experience a change in position angle but also in inclination. 
Thus if this H 1 1 region is located in the outskirts of the 
galaxy its position can be in the (warped) plane of the 
galaxy. 

The Ha distribution could, in theory, be severely 
affected by internal dust extinction, especially in the 
edge-on orientation. However, in the case of UGC 1281 this 
seems unlikely because no clear du st lane can be observed 
in HST imaging of UGC 1281 (|Bomans fc WeisI 120081 ') 
and dwarf galaxies ar e expected to have a low metallicity 
jPilvugin et al.l |2004| ) and therefore little dust content. 
The reddish color of UGC 1281 (|Makarova fc KarachentsevI 
1 19981 ) is most likely caused by the absence of a large young 
stellar population. This once more confirms the idea that 
UGC 1281 has a low SFR. 

In our data the maximum distance to the mid-plane 
where diffuse gas is still detected is ~25"(0.65 kpc). This 
is simila r to th e the extent of the stars (^ij25. mino r = 18", 
Ivan Zed l|2000l '). |d25, minor = 23", IdeV aucoule urs et al.l 
1 19921 ) ). ]f we fit an exponential to the inner vertical 
intensity profile of the ionized gas we find a scale height 
of 8.5" (0.22 ± 0.03 kpc) assuming that the galaxy is seen 
perfectly edge-on. When we follow the same procedure 
for the continuum emission in our spectra and an /-band 
imag43 we find a scale height of 7.6" (0.20 ± 0.01 kpc) and 
8.0" (0.21 ± 0.01 kpc) respectively. 



obtained through the NASA extragalactic database 



3.1.2 Hi Distribution 

Fig.[2]shows the DSS 2 red image of UGC 1281 overlaid with 
the contours of our integrated H I flux map. This H I map 
was constructed by adding all the channels of the Circular 
Beam data cube. To keep the addition of noise to a mini- 
mum, only regions that are above 3a in the Low Resolution 
cube (see Table [5} were considered. 

From Fig. [2] we can see that the H I of UGC 1281 is 
at first glance quite symmetrically and evenly distributed. 
However, a closer look reveals asymmetries and peculiar- 
ities in the H I distribution. It warps away from the ma- 
jor axis at about 90" (2.4 kpc) on the South West and at 
about 100" (2.6 kpc) on the North East. The warp initially 
shows the normal S-shap e observed in many edge-on galaxies 
(|Garcfa-Ruiz et al] |200!^l but bends back to the plane of the 
inner disk at larger radii. This behaviour is seen especially 
at the South West side, at a radial offset of ^200". T his 
warp was already observed by iGarcia^ Ruiz et all (|2002|) m 
the W HISP observations of this galaxy ( van der Hulst et al.l 
l200j ). 

When we compare the lowest contour in the integrated 
moment map of the WHISP observations with our own (see 
Fig. (2] red contour) , we see that in our observations more 
emission is detected in the radial as well as in the vertical 
direction. The growth in both directions is similar in extent 
and this indicates that even with our deep observations we 
might not be detecting the lowest levels of emission of this 
galaxy. 

Furthermore, the H I distribution displays a central de- 
pression. This depression appears as a ring around the center 
of the galaxy and ranges from 10"to 40" (0.26-1.05 kpc) ra- 
dial offset from the center of the galaxy. It appears to be 
symmetrical around the center of the galaxy. 

The H I in UGC 1281 shows significant extensions away 
from the major axis. The Hi extends up to 70"(1.8 kpc) 
on both sides of the plane at column densities Nhi = 
4.0 X 10^^cm~^(3a). This extent is much more than the 
FWHM of the beam (26") which is clearly seen in Fig. El 
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Figure 3. Intensity profiles of tlie observed Hi. Left figure: Ver- 
tical intensity profile, averaged over the inner 100"in the radial 
direction and normalized to the peak intensity. Right figure: Ra- 
dial intensity profile, averaged over the inner 40" of the galaxy 
in the vertical direction. Black line: Data, Green squares: Best fit 
model. Red diamonds (Blue triangles): best fit model with a scale 
height -2" (-1-4") (see §|4)l. Blue dotted dashed line: beam. 

This figure shows the vertical distribution of the data (black 
solid line) and a Gaussian (blue dot-dashed line) with a 
FWHM of 26". Both are normalized to the maximum of 
the data. In this figure it is easily observed that the wings 
of the data are much more extended than the observational 
beam. 

3.2 Velocity Distribution 

3.2.1 Hot Velocities 

Fig. U shows the velocity field of the PPAK observations. 
This velocity field was obtained by taking the peak position 
of the Gaussian fitted to each bin (see § 13.1. If) . This is by no 
means equal to the actual deprojected maximum rotational 
velocity in the galaxy but it is an apparent mean velocity 
determined by a combination of the rotational velocity, the 



Figure 5. Velocities along the major axis. The left panel 
shows the Ha velocities: Black symbols are the PPAK data pre- 
scnt ed in this paper, blue symbols are the data obtained by 
Kuzi o de Narav et al. I 1I2OO6I) . Right figure shows Hi velocities: 
Black symbols are the H I, the blue line shows the velocities ob- 
tained from the model. The red lines are the input unprojected 
rotation curves of the best fit models described in § |4] The arrows 
indicate the positions of the velocity cuts parallel to the minor 
axis in Fig. [G] 

density distribution of the ionized gas, and the opacity of 
the dust. From here on whenever we mention velocity we 
are referring to this mean velocity unless otherwise noted. 
The gaussian fitting procedure, and therefore the mean ve- 
locity, was chosen because with a channel separation of 70 
km s^^ it is impossible to confidently fit the intrinsic shape 
of the emission line. 

Fig. [S] (left) shows a cut 10"wide of the velocity field 
alo ng the major axis. Ove r plotte d are the velocities obtained 
by iKuzio de Narav et al.l l|2006l) with the DensePAK IFU 
(blue symbols) and the rotation curve ob tained from the 
modeling (see § |4]). iKuzio de Narav et al.l 12006) were not 
able to trace emission as far out in radius with their observa- 
tions. Since their exposure time and fiber size is equal to ours 
this is most likely caused by the fact that they do not bin the 
data and the lower sensitivity of the DensePak IFU. The ve- 
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Figure 4. Velocity field of the ionized gas. The field was constructed by taking the central position of the fitted Gaussian in all binned 
spectra. The systemic velocity (Vays = 156. km s~^) has been set to 0. The separate pixels show the fiber positions and the colors run 
from -60 to 60 km s"!. The arrows indicate the positions of the velocity cuts parallel to the minor axis in Fig. [S] 
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Figure 6. Four cuts parallel to the minor axis through the velocity field at -71, -23, 1, 25"radial offset from the center. The negative 
(positive) offsets are the So uth East (North West ) in the sky. Black points arc the unbinncd PPAK data. Blue, voronoi binned PPAK 
data. Red, Fll data. Green. TKuzIo de Narav et al. I f2006i ) data. 



locities obtained bv lKuzio de Narav et al.l (|2006l ) agree well 
with our values, which assures us that there are no system- 
atic errors in our reduction or the gaussian fitting procedure. 

In this plot we see clearly that the part of the galaxy ob- 
served in our Ha field of view is still resembling a slow rising 
rotation curve that indicates solid body rotation. This be- 
havior of the rotation curve is quite typical for dwarf galax- 
ies, which all seem to have a large inner region in which their 
rotat ion curve resembles solid body rotation l|Cote et al.l 
bOOOt ). This behavior could be caused by us considering the 
mean velocities. However, this is unlikely as we will discuss 
in §[3221 

Another thing that is quite clear is that the flux peaks 
in the Ha along the major axis stand out in the velocity 
curve as areas with lower velocities (Fig. [5l left, black ar- 
rows). This is expected if they are Hll regions, since the 
chance that they would lie on the line of nodes is low, and 
would imply that most of the peaks we see in the Ha distri- 
bution are indeed H 1 1 regions. Another possible explanation 
is that the intensity peaks correspond to a higher density in 
the radial density profile of the galaxy with thicker clumps 
offset from the line perpendicular to the line of sight. This 
explanation is supported by the fact that some of the peaks 
appear symmetrically around the center of the galaxy (see 
Section [3. l.ip and that the star formation rate in UGC 1281 



is so low (SFR = 0.0084 M0 yr."^) that giant Hll regions 
would not be expected to reside in the galaxy. 

If so, this would require that as the influence of the over- 
density or H 1 1 region, on the shape of the line profile dimin- 
ishes, the measured velocities should increase. We can see 
that this is the case when we look away from the mid-plane 
as is illustrated in Fig. [S] These plots show the velocity as 
function of vertical offset from the plane in bins 12"wide for 
the binned an d unbinned Ha (blue and bla ck, respectively), 
the data from lKuzio de Narav et al] l|2006l ) (green), and the 
H I (red) at the positions of the peaks|3 Here we see that 
the velocities move away from the systemic velocity as we 
look higher above the plane. This confirms the idea that the 
velocities on the major axis are lowered by a line-of-sight 
projection effect. These results and their implications will 
be discussed in Section [5.21 in combination with the results 
of the H I. 



^ Note that also above the major axis, the DensePak data and 
our data agree within the errors. 
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Figure 7. Velocity field of the neutral hydrogen. The field was constructed with the task MOMENTS in GIPSY (see text). Contours 
are from -60 km s~^ to 60 km with increasing steps of 10 km s~^. Black contours are the data. The white contours in the upper 
panel show the line-of-sight warp model. In the lower panel the white contours show the best fit edge-on model with a lag of 10.6 km 
kpc~l (see 



3.2.2 Hi Velocities 

Fig. [7] shows the velocity field of the H I observations. 
This velocity field was constructed with the MOMENTS 
routine in GIPSY by selecting the first moment of the data 
cube. The routine determines the intensity weighted mean 
velocity position of the peak of the line profile in every 
pixel of the cube. For symmetric profiles, this is analogous 
to fitting a gaussian profile to the emission line and taking 
its velocity at the peak. Since UGC 1281 has only small 
rotational values the line profiles are almost symmetric. We 
have checked this by comparing a map with the velocities 
where the line profiles have their maximum with this 
GIPSY map and we find no differences greater than 6 km 
s~^, which is less than the velocity resolution (see Table 
[T|. Thus, these velocities are in principle comparable one 
to one with the velocities of the Ha and differences should 
be due to the conditions of the gas (distribution effects, 
dust, resolution, self absorption, or a real difference in the 
rotational speed of the ionized and neutral gas). 

It is common to retrieve the rotation curve of edge-on 
galaxies through envelope tracing. Theoretically this is the 
correct way of retrieving the rotation curve of an edge-on 
galaxy (|Sancisi fc AllenI [1979I ). however there are several 
reasons why Gaussian fitting is preferable in an edge-on 
dwarf system. As explained in the previous paragraph the 
velocities in a dwarf galaxy are small, this means that the 
dispersion of the gas makes up a significant part of the line 
profile. Since in envelope tracing, either by fitting gaussians 
or a scaling of the maximum intensity l|Sofue fc RubinI 
I2OOII ). this dispersion is a chosen parameter it increases 
the uncertainty of envelope tracing. Furthermore, envelope 
tracing is also an estimation of the real rotational velocities 
and detailed modeling is still required. Therefore we prefer 
the well understood and described method of Gaussian 



fitting. This does mean that only from the modeling we 
obtain information about the real rotation curve. 

Fig. [5] (right) shows a cut 10"wide of the velocity field 
along the major axis. Here we see the same slow rise in the 
rotation curve as seen in the Ha but we also see that in 
the H I we do reach the flat part of the rotation curve at 
~120"radial offset, outside the Ha field of view, from the 
center of the galaxy at a maximum velocity ~60 km s~^. 

If we move away from the plane (Fig. [7] and [HJ we see 
that the velocities are lower than on the major axis. This 
implies that either the outer parts of the galaxy are heavily 
inclined, that the disk is flaring or that the gas above the 
plane is 'lagging'. We can quantify this vertical velocity 
gradie nt in normalized PV diagrams (see iKamphuis et al.l 
(|2007h . § 6.2) parallel to the minor axis (Fig.[8l bottom row) 
by fitting a straight line to the maxima between 20"and 
50"offset from the plane. Of course this fit is affected by the 
warp, lag, beam smearing and possibly other systematics. 
Therefore, the measured value has to be compared to the 
same measurements of models that can explain the observed 
declining velocities. These models will be presented in the 
next section and discussed in Section [S] 



4 MODELS 

The Ha observations show such an irregular distribution 
that constructing a model for comparison will not bring 
more insight into the structure of the Ha distribution. Also, 
the distribution is peaked in many places whereas our mod- 
els would resemble a smooth exponential distribution. We 
therefore decided to forego any attempt at modeling the 
ionized hydrogen in 3D. However, there are several parame- 
ters that can be derived from the data or obtained from the 
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Figure 8. Normalized PV diagrams parallel to the minor axis averaged over 100" at a radial offset of -175"to -75"(left) and 75"to 175". 
From top to bottom, the edge-on model without lag, the edge-on model with lag, the line-of-sight warp model and the data. The white 
contour outlines the 0.75 level of the data. The solid lines show the lag obtained from an edge-on lagging model (see ii 14.11 1. The dashed 
lines show the derived error on the lag. The zero point is set to the point of maximum intensity in a non-normalized PV diagram. 



Parameter 


Value 


Reference 


EM Scale height 


0.22 ± 0.03 kpc 


[11 


Scale length 


1.09 kpc 


[2] 


SFR 


0.0084 Mq yr. 


[2] 


Ha luminosity 


3.81±0.01 xlO"" erg s"! cm'^ 


[3] 



Table 4. Ho parameters that could be measured from t he obser- 
vation s or obtained from the literature. [1] this paper, [2] Ivan 
1I2OOII') . [3] Ivan Zed ll200d'l . 



literature. These are presented in Table [4] for comparison to 
the H I model. 



4.1 Tilted ring models 

To better understand the 3D distribution and the kinematics 
of the neutral hydrogen in UGC 1281 we have constructed 



a range of tilted ring models. These models were created 
with the GIPSY task GALMOD. The GALMOD routine 
was modified by one of us (G. H.) and F. Fraternali to be 
able to include radial and vertical motions as well as a ver- 
tical gradient for the rotational velocities. 

The modeling is not straightforward, since many differ- 
ent parameters may be independently fitted and some of 
them are degenerate. In the case of an edge-on galaxy mat- 
ters are even more complicated because small changes in 
density in the outer rings can seriously affect the inner rings 
of the model. For this reason the fitting always has to start 
at the outer radii and then move inward ring by ring. Since 
it is clear that the warp is not symmetric at all, we fit both 
sides of the galaxy independently. The goodness of a fit was 
determined by investigating several representations of the 
data (PV diagrams, residual maps) by eye. 

The following, iterative, procedure was followed for the 
fitting. We start by comparing radial and vertical density 
profiles of the data with the models (see Fig. [31). This has to 
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Figure 9. Upper panel: Integrated Hi flux map of UGC 1281 (see Fig.|2]l. The black contours are the data at the same levels as in Fig. 
[T] The red contours are the same levels for the best fit edge-on model. Lower panel: Residual map of the data minus the best fit edge-on 
model (total intensity). 



be done to get a initial guess of the density profile and scale 
height. When these profiles have a reasonable fit we start 
fitting the position angle on a integrated moment map. The 
position angle was deemed to fit when, at all radii, the ver- 
tical position of the peak value was equal in the models and 
the data. 

In the tilted ring model the inclination of each ring is 
related to the position angle and dependent on the angle 
between the axis on which the ring tilts and the line of sight. 
We assume a simple warp, where the tilt axii for all rings 
are aligned, e.g. a straight line of nodes. The assumption of 
a straight line of nodes is a simplifying one that lowers the 
degrees of freedom in the model, making it possible to con- 
strain the vertical distribution and kinematics in the model. 
By assuming a simple warp the inclination of each ring is 
coupled to its position angle by just one parameter, the an- 
gle between the line of sight and the line of nodes. We shall 
refer to this angle as the warp axis angle. 

From the data it is impossible to determine an exact 
warp axis angle. We therefore set out to find the maximum 
and minimum warp axis angle that can fit the observations. 
Starting at a minimum warp axis angle of 0°, e.g. all rings 
90° inclined, we found an acceptable density distribution 
immediately. 



We then set out to find the maximum warp axis angle. 
This maximum is set at the point where the vertical profile 
starts to show a break, which is not observed in the data 
(see Fig. |3} . By slowly increasing the warp axis angle and 
comparing the vertical distribution of the model to the data 
we find that this break becomes significant when the warp 
axis angle is greater than 60°, thus defining the maximum 
line-of-sight warp model to have a warp axis angle less than 
60°. 

At this stage we split our modeling into two parts. Be- 
sides a model with a warp axis angle of 0°, the edge-on 
model, we construct a second model with a warp axis an- 
gle of 55°, the maximum line-of-sight warp model. Now the 
only parameter of the density that is not determined yet 
is the vertical distribution of each ring. In our models the 
vertical density, of all the rings, declines as an exponential 
of which the steepness is set by its scale height. By increas- 
ing/decreasing the scale height of our models systematically 
we find the upper and lower limit of the scale height. We 
allow the scale height to change from ring to ring but only 
increase toward larger radii. This increase would resemble a 
fiare. In this way we find that the best fit scale height of the 
H I in UGC 1281 is a flaring model that ranges from 10"to 
15" (0.26 to 0.39 kpc) for our edge-on model and from 7"to 




Figure 10. Color plot of the Hi PV diagram along the major axis. Overlaid with contours of the data (grey/black) and the best fit 
model (red). The contour levels are -3cr, -1.5o", 1.5a-, 3(t and 6a {I a = 0.5 mjy/beam) etc. 




Figure 12. Left (Right): Color plot of the Hi PV diagram at 26" (-26") (±0.6 kpc) offset from the major axis. Contours are at -Scr, - 
1.5(T, 1.5 (T, 3 a, 6(T and so on. The black solid contour is the data, red contour is the best fit model with no lag, white contour is the 
best fit model with an assumed vertical gradient of 10.6 km kpc~^ in the rotation curve. The colormap ranges are indicated above 
the panel in Jy/beam. 



12" (0.18 to 0.31 kpc) for our maximum line-of-sight warp 
model. 

In order to estimate the accuracy of our scale heights 
we systematically raise and lower the best fit flaring model. 
It turns out that the vertical profile significantly deviates 
from the observations when we add (subtract) more then 
4" (2") (100 and 50 pc respectively) to the scale heights of 
our best fit models. The vertical distributions of these mod- 
els are shown in Fig. [3] as blue triangles (upper limit) and 
red diamonds (lower limit) for the edge-on model. The blue 
dot-dashed line indicates the spatial resolution (FWHM = 
26"). 



As a last check on the density distribution we construct 
an integrated moment map of the residual cube. Fig.[9]shows 
the residual map of the best fit density distribution. This fig- 
ure shows only the density distribution of the best fit edge-on 
model because, above the sensitivity limit of the data, the 
differences in the density distribution between the two mod- 
els are less than 5%. Therefore the difference between these 
two models would not be visible in such plots. In this figure 
the wiggles that necessitate adding a warp, purely in PA, 
to the edge-on model, and their reproduction in this model, 
are clearly seen. 

After we obtain a satisfactory density distribution we 
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Figure 13. Same as Fig. I12l but now at ± 56"(itl.3 kpc) offset from the plane. 



start fitting tiie rotation curve and dispersion of the gas. 
Initially we set the dispersion of the gas to be constant at 9 
km and then start fitting the rotation curve. As a first 
guess for a lower limit of the rotation curve we use the ve- 
locities from the Ha and H I. For the inner ± 50"we use the 
Ha velocities and beyond this point we use the velocities 
obtained from the Hi first moment map (see Fig. [5|. This 
rotation curve is then raised from the outside in until it fits 
the data. 

As an upper limit we start from a flat rotation curve and 
fit this curve by lowering the velocities from the inside out. 
The change between rings is always half of the difference 
between the previous two rings, with a minimal increase of 
4 km s~^ for the inner rings. Because a low inner density 
could artificially lower the observed velocities in an edge-on 
system we have set the densities of the inner rings to zero 
while fitting this rotation curve. 

Once more the fitting is an iterative process where the 
goodness of the fit is determined by comparing the ma- 
jor axis position-velocity (PV) diagram of the data to the 
model. An example of such an PV diagram is shown in 
Fig. IIOI where the best fit edge-on model is plotted with 
the data. The color scheme and the black contours are the 
data whereas the red contours display our best fit edge-on 
model. Again only the edge-on model is shown because also 
in this PV diagram the differences between the two models 
are less than 5%. When we have obtained satisfactory fits 
to this PV diagram for both the upper and lower limit rota- 
tion curves we find that they differ no more then a channel 
width (4.1 km s~^) at any ring position which indicates that 
the rotation curve is well constrained. Besides the upper and 
lower limit rotation curves matching very well, also the rota- 
tion curves of both, independently fitted, sides of the galaxy 
match up quite well. The final rotation curves are shown in 
Fig. [51 and 1141 (Top Panel). These curves are the average of 
the upper and lower limit of both best fit models. Fig. [10] 
show the rotation curve of the best fit edge-on model as a 
black dashed line overlaid on the major axis PV diagram. 

After fitting the rotation curve we found that the fit to 
the major axis PV diagram could be significantly improved 
by introducing a gradient in the velocity dispersion. This 
gradient runs from (Tv=11 km s~^ in the center to crv=8 km 
s~^ at the largest radii (see Fig. 1141 Bottom Panel). 

The model cubes should now be comparable to the 
data cube everywhere and any major deviations can only 



be caused by gas that is deviating from co-rotation at high 
projected distances from the mid-plane. In the case of the 
maximum line-of-sight warp model the observations are fully 
reproduced by the best fit model. However, in the case of the 
edge-on model the velocities above the plane are clearly over 
estimated. This can be seen in Figure [8] When we compare 
the edge-on model (Upper panels) to the data (Lower pan- 
els) it is easily seen that in the data the emission peaks at a 
vertical offset lie at lower velocities than the emission peak 
in the mid-plane. However, such a bending is only produced 
by the warp when the rings are not edge-on and thus com- 
pletely lacking in the edge-on model without a lag. Therefore 
we reproduce the edge-on model with a rotation curve that 
declines as a function of distance to the plane, e.g. a lag. 
We construct seven models where the vertical gradient is 
increased from 4 to 28 km s~^ kpc~^ in steps of 4 km s~^ 
kpc^V 

When we measure the vertical gradient in these models 
in the same way as in the data (see § I3.2.2p . we find that 
the value measured from the data corresponds to a lag of 
10.6 ± 3.7 km s~^ kpc~^ (see Fig. lllf) in the edge-on model. 
After adding this lag to the model it was found that, due to 
beam smearing effects, the major axis PV-diagram showed 
emission at slightly lower velocities than the data and the 
line-of-sight warp model. To correct for this effect the fiat 
part of the input rotation curve was raised by 2 km s~^. 
With this correction the major axis PV-diagrams of both 
models agreed within 5% again. A model with this lag pro- 
duces a satisfactory fit at all heights. 

As a last check to the model we compare the non-lagging 
model to the data and the lagging edge-on model in PV dia- 
grams at two distances above (below) the mid-plane. Figures 
[12] and [13] show PV diagrams parallel to the major axis at ± 
26" and ± 56"offset from the plane. Already Fig. [T^] is sug- 
gesting that a model with a lag gives a better fit, though we 
cannot exclude the non-lagging model at this height. At a 
height ~60" above the plane (Fig. I13p it is clearly seen that 
the velocities around the 3(t contours are too high in the 
non-lagging model. 

Figure [5] shows also the normalized PV diagrams paral- 
lel to the minor axis for the best fit lagging edge-on model. 
Also in these figures it is now seen that the edge-on model 
with a vertical gradient and the maximum line-of-sight warp 
model fit the data equally well. 

Fig. [14] shows the parameters for the best fit models 
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Figure 11. Input lag vs. Measured lag in the models. The solid 
line shows the fit to the measured values (displayed by crosses) of 
the models. The horizontal solid line displays the value obtained 
from the data. The vertical solid line shows where the measured 
value intersects the model fit. The grey shaded area indicates the 
error in the measurement and the dashed vertical lines where the 
error intersects the fit. This range determines the error bar on the 
value of the lag. 



of UGC 1281 (Edge-on model (black lines) & Line-of-sight 
warp model (symbols)). The best fit for the edge-on model 
has a flare and the scale height ranges from 10" (0.26 kpc) in 
the inner parts to a maximum of 15" (0.39 kpc) on the North 
Eastern outer parts. The change in PA is the same for both 
models but in the case of the line-of-sight warp model this is 
coupled to a change in inclination of the ring as previously 
explained. The scale length of this model galaxy is 46"(1.2 
kpc). The scale length is the same for the maximum warp 
model, which is also flaring and has a scale height of 7"(0.18 
kpc) in the inner parts and 12" (0.31 kpc) in its outer parts. 



4.2 Ballistic Models 

To investigate what we would expect for UGC 1281 in the 
sense of extra-planar gas brought up from th e disk by su- 

Jiernovae we have used the ballistic model of ICoUins et al.l 
2002[ ). In this model gas is blown out of the disk into the 
halo, with an initial vertical velocity (called the kick ve- 
locity). The model naturally predicts a vertical gradient in 
rotational velocity which has a higher magnitude in model 
with a high kick velocity. To obtain the vertica l gradient of 
the m odel we follow the procedure outlined bv lHeald et al.l 
l|2007h , with the following difference. The initial disk for the 
model is infinitely thin, this is because the scale height of the 
Hq in UGC 1281 is already comparable to the scale height 
normally used for the initial disk (0.2 kpc). Therefore this 
model should be thought of as an absolute upper limit on 
the kick velocity. 

We find that the scale height of 0.2 kpc is reproduced 
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Figure 14. Model parameters as function of radius. For the edge- 
on model (Solid line: South West side. Dashed line: North East 
side) and the maximum line-of-sight warp model (Blue plusses: 
South West side. Red crosses: North East side). 



when the initial kick velocity is 10 km s~ . We run the sim- 
ulation until the system reaches steady state (~1 Gyr). At 
this point there are no clouds in the halo at radii smaller 
than 4 kpc due to the radial outward movement of the 
clouds. Therefore we cannot measure the vertical gradient in 
the rotational speed at these radii. However, at radii between 
5-10 kpc we find a small gradient of 0.4 km s~^ kpc~^. This 
is a very low value as one would expect for a small galaxy 
like UGC 1281 since kick velocities must be small in a low 
potential, because otherwise the scale height becomes too 
large. 
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5 DISCUSSION 

5.1 Hydrogen in the plane 

When we compare the parameters of the best fit H I models 
to the parameters extracted from the Ha data (Table [4]) 
we see that compared to the ionized gas the neutral gas is 
much more extended in the radial direction. Apart from 
this its scale length is slightly larger. So, the ionized gas is 
more centrally concentrated than the H I. 

Note that both the neutral gas and the ionized gas 
display a central depression, or peaks symmetrically around 
the center, in their distribution (see Section 13. If) . Also, 
the H[q distribution is more irregular than the neutral 
hydrogen. This can easily be seen in Figures [1] and [2] by 
counting peaks in the distribution (see Section [3.1.ip . One 
of these 'peaks' in Ha has a clear offset from the plane of 
the galaxy and this region does not stand out in the H I 
observations. 

Another effect is the start of the warp. If we look 
at the integrated moment maps the neutral gas starts to 
deviate from the major axis at a radius almost two times 
larger than the radius where the Ha bends away from the 
major axis (50"and 90", respectively). Most likely this is 
a resolution effect, because when we look at the best fit 
model for the H l we see both the H I and the Ha start to 
warp around 50" (1.3 kpc). 

When we compare the measured velocities of the 
neutral and ionized gas they look very similar on the major 
axis (see Fig. [5]). Ideally this should only occur if the 
sizes of both disks are similar and there is no absorption, 
self-shielding or dumpiness affecting the emission. We know 
from the analysis of the Ha that the velocities in the plane 
are lower than those above the plane in several places. Also, 
the H I velocities are slightly lower, due to beam smearing, 
when compared to the best fit model. Since for the Ha 
emission we can only look at the rising part of the rotation 
curve, where the rotational velocities are low, deviations 
from the real rotational velocities are small. Therefore, 
it is not unreasonable that both effects provide us with 
the same velocities since all absolute deviations from the 
rotational velocities are small. Fig. [5] (Red line) shows us 
that the real rotation curve, obtained from the modeling, 
is slowly rising and reaches a flat part around 60 km s^^. 
This means that UGC 1281 has the slow rising rotation 
curve which is common in dwarf galaxies, but it also shows 
clear differential rotation. 



5.2 Hydrogen above the plane 

We see (Fig. [T]) that the Ha extends less far in the vertical 
direction than the H I. Despite the fact that a difference in 
sensitivity may cause us to miss Ha emission at distances 
further from the plane than 20", the scale height, measured 
by fitting a single exponential to the vertical distribution, 
of the Ha is, equal at best, but most likely smaller than the 
H I scale height. 

Even though the vertical extent of the Ha is limited 
and not supplying much information above the plane the 
H I extends well above the plane and much can be learned 
from it. If we assume a constant gas dispersion of 9 km s~^ 



we would expect a scale height of ~400 p c for the disk, 
from a simple estimate (jPuche et al.l 1 19921 ) . This implies 
that the gas at large distances from the mid-plane is not an 
additional component but merely the high-latitude gas of a 
pressure supported disk. This idea is strengthened by the 
fact that a ballistic model requires a very low kick-velocity 
(~10 km s~^), of the order of the velocity dispersion in the 
model, and indicates that the gas at large distances from the 
mid-plane is not pushed to high latitude by star formation 
processes. Furthermore, continuum maps constructed from 
the line- free channels of the 21 cm observations do not 
show any continuum emission at the position of UGC 1281 
once more confirming the low SFR. This implies that UGC 
1281 does not contain a typical halo such as seen in NGC 
891. However, also in the case of UGC 1281 the kinematics 
show that as the distance from the mid-plane increases the 
projected mean velocities decline (see Fig. [8|. 

From our modeling we obtain two possible explanations 
for the high latitude H I and its declining projected mean 
velocities. We find that a maximum line-of-sight warp 
model, where we push the warp as much into the line of 
sight as the data allows, fits the data equally well as a 
purely edge-on model with a warp in the plane of the sky 
and a vertical gradient in its rotation curve of -10.6 ± 3.7 
km s~^ kpc~^. Both models assume a single exponential in 
the vertical direction and the edge-on model does assume 
that the radial gas distribution above the plane is similar 
to that in the disk. In theory the lower velocities could 
also be caused by an alternate distribution above the 
disk. However, since the vertical gradient is seen almost 
everywhere in the observations, it seems unlikely that this 
is the case unless the gas is in a warp or flare and these 
possibilities are included in the modeling. 

To outline the importance of separating between a lag 
and a line-of-sight warp we will discuss these two options 
and their implications in two separate sections below. These 
discussions are by no means meant as a complete overview 
of the theory and observations behind extra-planar gas but 
only to highlight the parts where UGC 1281 can signifi- 
cantly contribute to a better understanding of the theories. 
For a full and comple te review of cold e xtra-planar gas 
we refer the reader to ISancisi et al.l ()2008l ) and references 
therein. 



5.2.1 A warp 

The first possibility that can explain the observations is 
that the gas at large projected distances from the mid-plane 
in UGC 1281 is located in a simple warp. Warps are quite 
common in disk galaxie s and often they are asymmetric 
(|Garcia-Ruiz et ahlliooi ) as is observed in UGC 1281. If 
the H I and its kinematics in UGC 1281 are to be fully 
explained by a warp, the warping axis has to have an angle 
of 55° with respect to the line of sight. 

We would like to point out the very small differences 
between the edge-on model and the maximum line-of-sight 
warp model. The max;imum difference in inclination between 
the two models is 10° but on average 5° in inclination. This 
results in a difference in scale height of ~3"(see Fig. I14p . 
These differences take away the need for a vertical gradient 
completely. This once more shows that differences between 
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line-of-sight warps and lagging halos are very subtle, and 
that great care must be taken to exclude one of the two 
options. 

When we measure once more the maxima in normalized 
PV diagrams parallel to the minor axis we find an apparent 
vertical gradient of 4.9 ± 1.1 km s~^ kpc~^, for the 
maximum line-of-sight warp model, which is consistent with 
the measurement from the data (6.9 ± 2.0 km s~^ kpc~^). 
Therefore there is no need to introduce a lag or other extra 
kinematical effects into this model. 

The best fit model for this case has central scale heights 
which are similar to those measured from the stars and 
the Ha. This would mean that the stars and the ionized 
gas hardly extend into the warped outer regions of the 
disk and that all the H I at large projected distances from 
the mid-plane is in the warp and fiare. We have already 
seen that the PA starts changing well within the maximum 
radius of the Ha as well as inside the optical radius. This 
is seen for every model, whether it is edge-on or maximum 
line-of-sight. 

The start of the warp could be affected by a slight error 
in the assumed PA. We have tested this by rotating the 
integrated H I moment map and the red DSS image (See 
Fig. 0) by an additional 1 and 2 degrees (PA= 39°, 38° 
respectively) and plotting the vertical offsets of the peak 
of the vertical profiles. This shows that for the red optical 
image the PA of 40° gives the flattest central distribution 
but that the start of the warp in the H I can be pushed 
outward by ~ 15" by assuming a PA of 39°. This means 
that the observed start of the warp remains well within 
the optical radius (iD25 = 134") even if the assumed PA is 
slightly ofi'. 

The fact that the warp starts wit hin the optical radius 
is inconsistent with the findings of Ivan der KruitI \2mi ) 
that warps start just beyond the truncation radius of the 
stellar disk. Even more so, we know from the observed 
Ha emission that also the ionized gas is slightly warped 
(see § I3.1.ip . This would imply that if the warp in UGC 
1281 is formed by the accretion of gas from the IGM its 
initi al disk was not rigid enough to stabilize the infalling 
gas (Ivan der KruitI [20071 ). If so, one would expect a clear 
difference between the onset of the warp, with respect to 
its truncation radius, between dwarf galaxies and massive 
galaxies. 



5.2.2 A lag 

In the case of a warp which is purely in the plane of the 
sky the modeling indicates a larger scale height of the H I 
than the scale heights measured from the stars and the 
ionized hydrogen. This vertical density distribution can be 
modeled with a single exponential, however the kinematics 
indicate the need for a vertical gradient in the rotation 
curve when the warp is purely perpendicular to the line of 
sight. As in the case of the super-thin galaxy UGC 7321 
l|Matthews fc Wood|[2003l ) the origin of this lagging neutral 
hydrogen gas in a Low Surface Brightness (LSB) dwarf 
galaxy, with low star formation rates, would be puzzling. 

If the high latitude gas in UGC 1281 is lagging this 
would have some implications to current theory. iHeald et"al] 
l|2007h have compared the lag, or the vertical gradient in the 



rotation curve in three massive galaxies. They find tentative 
evidence that when they scale the lag with the observed 
Ha scale heighl[f| this new parameter (dV/dhz) is roughly 
constant at about 20 km s~^ scale height"^. However, 
the galaxies compared are of similar mass. The results 
presented in this paper give us now the opportunity to take 
a tentative look at a class of galaxies with much lower mass. 
The dynamical mass of UGC 1281 is 6.3x10^Mq, measured 
at our last point of the rotation curve, as opposed to 
^ 1 X 10^^ M0 for the galaxies in the study bv lHeald et al.l 
i2003). 

Within the picture of a lagging disk, we would have 
dV/dhz = 4.7 ± 1.7 km s~^ per scale height. This value 
would be inconsistent with a constant dV/dh^ ~20 km 
s^^ which was found for the 3 massive galaxies. As shown 



by the ballistic models (i) I4.2|l . the lag expected purely on 
gravitational grounds would be much shallower than the 
one than the one that we include in this model. Such a lag 
would thus requir e addition a l effec t to be at play s uch a s 
described by e.g. iBeiiiaminI l|2002l ): iBarnabe et all l|200d ): 
iFraternali fc BinnevI (|2008| ). In any case, star formation in 
the disk would be of negligible influence on the vertical 
gradient. 

The previous discussion clearly shows the need to 
unambiguously determine a lag in a small galaxy such as 
UGC1281 as well as the need for larger sample of galaxies 
with quantified vertical velocity gradients. In the present 
case, due to the fact that the addition of a lag to the model 
of UGC 1281 does not significantly improve the match with 
the data, we prefer the conceptually simpler line-of-sight 
warp model. 



6 SUMMARY 

We presented 21 cm and Ha emission in the edge-on dwarf 
galaxy UGC 1281. This is the first time such sensitive H I 
data have been presented for a dwarf edge-on. 

The integrated Ha velocity map (Fig. [T| shows a non- 
smooth distribution on the major axis with several peaks. 
One of these peaks is actually located beneath the major 
axis. It is unclear whether this Hll region is located above 
the plane of the galaxy or in its warped outer parts. 

The integrated H I velocity map (Fig. \^ shows a quite 
regular distribution with a central depression. This central 
depression appears to be symmetrical in position around 
the center but from modeling it follows that it is somewhat 
deeper on the NE side of the galaxy. Such a central depres- 
sion is not uncommon for dwarf galaxies. 

Furthermore this map shows that UGC 1281 is warped 
in its outer parts and this warp resembles a "normal" S- 
shape at its start. However, at large radii the warp bends 
back towards the inner plane of the galaxy. 

For the interpretation of the kinematics of the high lat- 
itude H I gas we constructed velocity maps from the Ha 
and H I data. Also 3-D models with a modified version of 
GALMOD are constructed in GIPSY. This modified version 



^ Note that lHeald et al" I ll2007t) use electron scale heights instead 
of the emission measure scale heights. 
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enables us to construct kinematic models with a vertical gra- 
dient. 

The velocities obtained from the data show a slow rise 
in the inner part. This is also seen in the rotation curve 
obtained from the modeling and therefore unlikely to be an 
effect of the H l distribution or resolution. At about 120" the 
rotation curve flattens off to a maximum rotational velocity 
~ 60 km s~^. This slow rise is common for dwarf galaxies. 

From our modeling we find that our data are not sen- 
sitive enough to distinguish between a lag or a line-of-sight 
warp. Both models fit the data equally well and there is 
only a small difference between the input parameters of the 
models. However, the models do start to deviate at emis- 
sion levels slightly lower than our current sensitivity limit. 
Therefore one would most likely be able to separate between 
the two models with deeper observations. 

In the case of lagging high latitude gas the low verti- 
cal extent and the low fiux level of the Ha emission would 
indicate that this high latitude H I does not originate from 
galactic fountains. The H l scale height implies a normal 
pressure supported disk and thus no need beyond turbu- 
lence for a mechanism to bring it up from the mid-plane. 
If in such a disk the pressure, of the gas, is solely depen- 
dent on the density (e.g. barytropic) one would expect it 
to be co-rotating. However, we find in our analysis that the 
high latitude gas in this case has a lag of 10.6 ± 3.7 km s~^ 
kpc~^ when compared to tilted ring models. This lag could 
be caused by infalling gas or pressure gradients above the 
mid-plane. 

In the case of a line-of-sight warp the ionized hydrogen 
and the distribution of the stars would, for the most part, 
not extend into the warped region of the disk. However, the 
scale height in the central parts would be the same for the 
stars. Ha and H I. 

Regardless of which model fits the data best, maximum 
line-of-sight warp or edge- on, the warp starts well with the 
optical radius (D^.'i=4.46'. Ide Vaucouleurs et al.l l|l992h '). at 
a radius of ^ 50", which is unlike more massive galaxies 
l|van der Kruitll2007l ). 

The small differences in input parameters between a 
model with a lag and one with a line-of sight warp show 
that great care must be taken to distinguish between lag- 
ging halos and line-of-sight warps since small changes in the 
modeling can have a great effect on the velocity field. 
Summarising our main conclusions: 



• The rotation curve of UGC 1281 is slowly rising in its 
inner parts and flattens off to a maximum rotational velocity 
-60 km s"^ at 120" (3.14 kpc). 

• The neutral hydrogen in UGC 1281 is more extended 
than the stars and the ionized emission in its radial distri- 
bution. 

• The gaseous warp start well within the optical radius. 

• Our observations can be fitted by both a vertical gra- 
dient in the rotation curve and a line-of-sight warp. The 
observations are not sensitive enough to separate between 
these two options. However, the line-of-sight warp model is 
conceptually simpler and therefore preferable. 
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